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ABSTRACT 

In this report, models are developed for the apparent 
temperature and backscatter coefficient of vegetated ter- 
rain to illustrate the effects of vegetation on the sen- 
sitivity of these parameters to variations of soil mois- 
ture. Three types of terrain are simulated for both the 
passive and the active case: a uniform canopy over a 

smooth surface, plant rows on a smooth surface, and plant 
rows on a rough surface. In each case the canopy is de- 
fined by its overall dimensions and by its electric per- 
mittivity, which is determined from the Weiner model for 
dielectric mixture. Emission and scattering from both the 
soil and the canopy are considered, but atmospheric effects 
are neglected. 

The expression for the apparent temperature of vege- 
tated terrain includes terms for the apparent temperature 
of the bare soil and of the canopy. The bare soil term is 
modified by an exponential term to account for attenuation 
of radiation from the soil by the canopy. The term which 
represents radiation by the canopy assumes that the canopy 
is homogeneous layer. The expression for the apparent 
temperature of row crops is the average of the apparent 
temperatures of covered or shadowed soil and visible bare 
soil. A rough surface is represented by a surface rough- 



Ill 


ness factor and a coefficient of effective area of specular 
refraction. 

The expression for the. backscatter coefficient of veg- 
etated terrain is similar to the expression for apparent 
temperature. The term for backscatter from smooth, bare 
soil is determined by the physical optics method. The 
term for backscatter from the canopy is a modified form of 
a model for scattering from long, thin dielectric cylinders. 

Caluclated data indicate that the sensitivity of the 
apparent temperature and backscatter coefficient to varia- 
tions of soil moisture, decreases as the amount of vege- 
tation increases. It is shown that the same effect results 
from increasing signal frequency or angle of incidence. 

Several sets of equivalent vegetation states, defined 
as different combinations of parameters which yield equiv- 
alent data, are tabulated. It is shown that in many cases 
the height-density product of a uniform canopy may be used 
to identify equivalent states. 
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CHAPTER I 
INTRODUCTION 

The Potential for Microwave Sensors 
in Agricultural Remote Sensing 

Remote sensing techniques have become established 
tools for agricultural planning and management. However, 
only in recent years has the potential value of microwave 
sensors in agricultural applications been examined. Radar 
and microwave radiometers have the potential for perform- 
ing certain tasks presently performed by photographic sen- 
sors, and other tasks which cannot be accomplished through 
photography. For example, since the development of syn- 
thetic aperture radar systems, it has been possible to 
produce radar imagery with resolution of several meters 
from high altitudes. Such imagery may provide as much in- 
formation about certain types of terrain as photographic 
imagery. Furthermore,, the ability to obtain radar imagery 
does not depend on the time of day and is almost indepen- 
dent of weather conditions. 

Microwave sensors are valuable not so much for the 
characteristics of the instruments, but for the region of 



the electromagnetic spectrum in which they operate. The 
microwave band extends from 0.3 GHz to 300 GHz. However, 
most research in agricultural remote sensing has been con- 
fined to frequencies less than 36 GHz (centimeter wave- 
lengths] because of economic and technical advantages [1]. 
Radiation at the lower frequencies can penetrate clouds 
and fog without serious degradation of the signal. Using 
cross-polarization techniques, radar echoes from the ground 
can be detected even during rain. Furthermore, at the 
lower frequencies the depth of penetration through vegeta- 
tion and soil may be several centimeters. A discussion of 
microwave penetration of soil is given in Appendix A. 

The use of microwave frequencies for remote sensing 
does have certain inherent disadvantages. Data are often 
difficult to interpret because several variables affect 
the measurements: frequency, polarization, angle, permit- 

tivity and roughness of the surface, etc. Data from vege- 
tated terrain depend in general upon the characteristics 
of both the plants and the soil. Sometimes the total ef- 
fect of the unknown surface parameters is unpredictable, 
in which case no information about the surface can be de- 
rived from the data. However, if a single parameter has an 
overriding effect on the data, then variations of this 
parameter can be detected. 

Experimental proof of the variation of the 
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permittivity o£ soil with moisture content has stimulated 
research in the area of measuring soil moisture with micros 
wave remote sensors. Some radiometer data from bare 
ground indicate that the moisture content of the soil can 
be determined with fair accuracy from the apparent tempera- 
ture of the soil. However, similar data from vegetated 
fields show very little correlation with apparent tempera- 
ture. 

Obj ective 

There is no question that the permittivity of soil is 
a function of moisture content. Therefore, the moisture, 
content affects the emission and scattering of microwave 
energy from the ground to the same degree to which the per- 
mittivity affects these processes. Obviously, as the 
height or density of vegetation increases, the permittivity 
of the soil becomes less significant. The question is, can 
the moisture content of a soil surface which is more than 
slightly vegetated be measured with microwave sensors? 

The objective of this thesis is to determine the ef- 
fect of vegetation on the ability to measure moisture con- 
tent by microwave remote sensing techniques. Both active 
and passive sensors are considered. Apparent temperatures 
and scattering coefficients are computed as functions of 
soil moisture for several simulated vegetated terrains. 
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Scope of Report 

The models developed for calculation of the apparent 
temperature and scattering coefficient are intended to be 
only first approximations of vegetated terrain. To account 
for attenuation by and radiation from the vegetation, the 
canopy is assumed to be a homogeneous medium. Scattering 
from the canopy is accounted for by using a model for scat- 
tering from a particular class of vegetation. 

Since the investigation is concerned only with the 
effects of vegetation and soil moisture on emission and 
scattering, the effects of atmospheric propagation have 
been ignored. A discussion of atmospheric effects is given 
in Appendix B. 

Chapter II briefly reviews some theories and experi- 
ments which provide the background for this investigation. 

In Chapter III models are developed for the apparent 
temperature of some simple vegetated surfaces. To include 
the effects of vegetation, the permittivity of the plant 
canopy is derived from a model for the permittivity of a 
dielectric mixture. Expressions are derived for the ap- 
parent temperature of a smooth, uniformly vegetated surface 
and for row crops on smooth and rough surfaces . 

Chapter IV treats scattering from the surfaces modeled 
in Chapter III. Scattering from an arbitrarily rough sur- 
face is described by defining an effective area of 
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specular reflection. 

Chapter V presents microwave data obtained from vege- 
tated terrain, and compares them with the results of Chap- 


ters III and IV. 
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CHAPTER II 

THEORETICAL AND EXPERIMENTAL BACKGROUND 

The scattering and emission of electromagnetic waves 
from natural surfaces are complicated processes which are 
very difficult to describe because the characteristics of 
a natural surface are not easily quantified. Vegetation is 
even more difficult to describe than soil because its con- 
figuration is constantly changing. 

Despite the difficulties encountered, progress has 
been made in describing the parameters which effect scat- 
tering and emission from natural terrain and in developing 
practical applications of the theory. In particular, much 
effort has been devoted to developing techniques of re- 
motely monitoring soil moisture content. This chapter 
gives a brief review of some of the theories and experi- 
ments which have contributed to an understanding of the 
problem. 


Peake’s Model for Apparent Temperature 

The model for apparent surface temperature developed 
by Peiike [2] expresses the emissivity of a surface in terms 
of its scattering coefficients. By employing an empiri- 
cally derived expression for the differential scattering 
coefficient, the problem of determining the emissivity of a 
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surface for all angles is reduced to determining a measur- 
able constant at any particular angle. Peake's theoretical 
predictions are in fair agreement with experimental mea- 
surements of the apparent temperature of asphalt and grass. 
A brief description of Peake's development is given below. 

For an element of surface area S upon which radiation 
of intensity I© is incident from a direction ( 0 o, 4 >o), the 
differential scattering coefficient is. defined by 


y(l9* .^ 4 ,) = I, 

isc^a 


(II-l) 


where is the intensity of radiation at a distance R from 

the surface in the direction (0 * d) ) . 

s' 

Considering that either of two orthogonal polarization 
states may be specified for both the incident and scattered 
radiation, there are actually four scattering coefficients. 
Each of these may be represented in general . 

The subscript i indicates the polarization of the incident 
radiation, and j the polarization of the scattered radia- 
tion. The letters o and s represent the directions (0©, 

00 ) and ( 9 ^, 0 ^). As a consequence of the reciprocity 
theorem, the scattering coefficients satisfy the relation, 




( 11 - 2 ) 
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where each of i and j may represent either polarization. 

The albedo and the absorption coefficient of the sur- 
face are both determined by the differential scattering 
coefficients. The albedo is defined as the fraction of 
power (of a specific polarization) incident on a surface 
from a direction (0o, «J>o) which is scattered into the upper 
hemisphere. Considering that the scattered radiation may 
be of either polarization, the albedo is 




The absorption coefficient a^(o) is the fraction of power 
incident from direction (6o> <po) with polarization i, which 
is absorbed by the surface. Considering only natural sur- 
faces which are so thick that no transmission occurs, the 
absorption coefficient is simply 


ajo) = I - Aj(o) 


(II-4) 


The emissivity of a surface can be related to the 
absorption coefficient by deriving a generalized statement 
of Kirchhoff's law. The emissivity is defined as 
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e- (o) 

Considering a surface in thermodynamic equilibrium with the 
upper hemisphere, it is assumed that as much energy of a 
given polarization leaves the surface in a given direction 
as falls upon it from that direction. This assumption 
leads to 



Applying the reciprocity relation of (II-2), (H-5) be- 
comes 

I = ejo) + +y..{o,s)'jdD.^ cii-6) 

The last term of the expression is now recognized as the 
albedo of the surface. Comparing (II-6) with [II-4], it is 
clear that the emission coefficient, for either polariza- 
tion, is equal to the absorption coefficient. 

The surface emissivity determines only one of three 


Power emitted with polarization i 
by unit area of surface into element 
of solid angle dJ^o in the direction 6o4o 

Power emitted with same polarization by 
unit area of blackbody at same tempera- 
ture into same element of solid angle 
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components of the apparent temperature. The total apparent 
temperature for a particular polarization may be expressed 
generally as 



T + T + T 
g s p 


(II-7) 


where = contribution of the thermal radiation emitted 
by the surface 

Tg = contribution of scattered diffuse radiation of 
the atmosphere 

Tp = contribution of scattered radiation from point 
or quasi-point sources 

A complete model of the measured apparent temperature 
should account for the effects of the atmosphere between 
the ground and the sensor. However, Peake's model assumes 
that the radiated power is measured near enough to the 
ground so that effects of the intervening atmosphere are 
negligible. 

From the definition of the- emission coefficient, the 
first contribution is simply the product of the emission 
coefficient and the actual ground temperature. 



Tg(o) = eCo)Tg 


(II-8) 


1 1 


The second contribution, that of the scattered diffuse 
radiation of the atmosphere, depends upon the scattering 
coefficients of the surface and the apparent temperature of 
the radiation incident on the surface from above. Assuming 
that T^(s) and Tj (s) are the apparent temperatures of the 
radiation with polarization states i and j incident from 
the direction (0 , <p ) , the contribution to the apparent 
surface temperature is 



(II-9) 


From the reciprocity relations it follows that 


TsiW - + T^(s)7ij(o.s)]o<nj cn-io) 


The third contribution is due primarily to the sun. 

Peake showed that for band the value of T . is much less 

u pi 

than one degree Kelvin. Therefore, the last term of (11-7) 
can generally be neglected. 

Writing the emission coefficient in terms of the 
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scattering coefficients, and assuming that = Tj(s), 

the apparent surface temperature may be expressed as 

Evaluation of (II-ll) requires that the scattering 
coefficients be expressed analytically. Peake found the 
following expression to be applicable and consistent with 
measurements . 




C^9o 

2c^e, 


( 11 - 12 ) 


where yo is a constant, v denotes vertical polarization, 
and h denotes horizontal polarization. The value of yo is 
on the order of 0.1. 

The advantages of Peake's model are that it is appli- 
cable to any surface satisfying the assumption of thermal 
equilibrium, and it is relatively simple to evaluate. No 
distinction is made between bare and vegetated surfaces. 
Each is considered simply a natural surface with a charac- 
teristic yo. 

The assumption of thermal equilibrium is generally 
violated by vegetated surfaces. Kumar [3] has reported 
that the temperature of plants may be as much as 20°C above 
ambient. Furthermore, plants may warm or cool within 
minutes, or even seconds. Nevertheless, Peake's model has 
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been used to predict with fair accuracy the apparent tem- 
perature of certain vegetated surfaces [2], 

Peake's Model for Long, Thin Dielectric Cylinders 

To better predict the scattering from certain types 
of vegetation, Peake developed a model for scattering from 
long, thin dielectric cylinders [4]. The cylinders are 
assumed to be randomly spaced and oriented, with greater 
probability of vertical orientation. Since the model 
accounts for attenuation and multiple scattering within the 
plant canopy, it seems to predict the scattering coeffi- 
cient of many vegetated surfaces rather well. However, 
application of the model is limited by several . assumptions 
or restrictions, as described below. 

All of the cylinders are assumed to be semi-infinite 
in length, extending downward from the plane z = 0. Con- 
sequently, scattering from the ground need not be consider- 
ed. 

All cylinders are assumed to be identical in shape, 
size, and electrical properties. The cylinders are round 
with diameter much smaller than a wavelength. 

The incident field entering the volume which contains 
the cylinders is assumed to be attenuated exponentially as 
it propagates vertically into the medium. This attenuation 
is due to both absorption and scattering. The attenuation 
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constant is denoted a. 

The cylinders are randomly spaced and oriented. It is 
assumed that all azimuth angles are equally likely. The 
probability that a cylinder is directed at an angle 0^ from 
the vertical is p(6^)dn, where dn is a unit solid angle, A 
probability density 

pl^i) = c^^$- CII-13) 

is assumed. The average number of cylinders per unit area 
is denoted N. The spacing between cylinders is considered 
large enough with respect to a wavelength that the phase 
difference between scattered fields is random. Conse- 
quently, the total average scattering cross section is the 
sum of the individual cross sections. 

With the above assumptions, Peake derived four expres- 
sions for the bistatic cross section per unit area of vege- 
tated ground. Since these expressions are rather long, 
only the simplest is repeated here as an example: 





(11-14) 
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where 


Bp = 



|AN)(Ak*)[l 


1 + 1 

1^/1 

|2* 


f 28TT 

' 3 1 
. 5* ’ 

t a] 

(kJ 

+ 3(c0«2^i • 




2 Q-^ aA/v\Q^ Cd^ 4^s 


> 




In these equations, A is the cross sectional area of a cyl 
inder, N is the average number of cylinders per unit area, 
and + je^'* is the relative permittivity of each cyl- 
inder. 

Peake has estimated the permittivity of vegetation as 

€r' + J€/' = 2.5fl -f) + (11-16) 

where f is the fraction of water by weight in the plant, 
and is the relative permittivity of water. For micro- 
wave frequencies the latter parameter can be approximated 
by 


75 




(11-17) 


The factor t appearing in (11-14) is the transmission 
coefficient for the surface of a cylinder at normal inci- 
dence . It is taken to be 


t - 




(11-18) 
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The value of the attenuation constant appearing In 
depends upon the polarization of the incident radiation. 
This parameter is not easily evaluated in practice. How- 
ever, Peake has estimated the value of a for horizontal 
and vertical polarization as follows: 


a 


h 





( 11 - 19 ) 

‘ 1 ] 


As with all scattering models, there are certain ad- 
vantages and disadvantages of applying Peake's model. One 
of the advantages of the model is that it accounts for most 


of the factors known to effect the scattering cross section 
of vegetation; wavelength, incidence and scattering an- 
gles, water content, etc. However, the fact that the model 
accounts for all of these presents a disadvantage in that 
some of these parameters are difficult to evaluate accu- 
rately. 


A serious limitation on the application of this model 
is the assumption that the power scattered from the ground 
is entirely negligible. The model is applicable only to 
cases for which the height and density of vegetation are 
greater than certain minimum values. Peake's model cannot 
be applied directly to the present investigation, in which 
the variation of scattering cross section due to variations 
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of soil moisture content is of interest. 

Measurement of Permittivity of Soil 

Since the variation of permittivity with moisture con- 
tent is responsible for changes of the scattering cross 
section and radiometi*ic temperature of soil^ the behavior 
of the permittivity with varying moisture is of interest. 
Several investigators have measured the relative permittiv- 
ity of soil at different microwave frequencies. The pro- 
cedures used and results obtained by two of , these are dis- 
cussed below. 

Measurements by Lundien 

Measurements were made by Lundien [5] to determine 
the effects of frequency, density, moisture content, and 
soil type on the dielectric constant of soil. Relative di- 
electric constants of more than 400 samples of twelve types 
of soil were determined. The moisture content of the sam- 
ples varied from 01 to 51.5?;. Measurements were made at 
frequencies of 1.074, 1.125, 1.311, 1.412, and 1.499 GHz. 

The technique used by Lundien employs an L-band inter- 
ferometer system as illustrated in Figure II-l. Prior to 
making a measurement, a reference standard is placed be- 
tween the antennas, and the attenuator and phase shifter of 
the known leg are adjusted to give a null. The reference 




Figure II-l. Block diagram of L-band interferometer [5]. 

I—* 

00 
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standard is then replaced by a sample of known thickness, 
and the circuit is again adjusted to give a null. The per- 
mittivity of the sample can be determined from the differ- 
ence in attenuation and phase of the sample and the re- 
ference standard. 

Inaccuracies due to multiple reflections and edge ef- 
fects can be reduced or eliminated by placing the sample at 
an oblique angle with respect to the propagation path. 
Lundien used 4S°, but the selection of angle is somewhat 
arbitrary. 

Lundien determined the relative dielectric constant 
and the loss tangent of each sample from tables computed 
from the following relation: 



kxp( 

~S^) 

\-r‘ 

‘ exp 

H 24>] 


( 11 - 20 ) 


where T = complex voltage transmission coefficient 
T' = transmission phase 
r = reflection coefficient 


For perpendicular polarization, 

ca<t.d - 
^ O0<l9 + 

and for parallel polarization. 


(Il-Zla) 


CII-21b) 
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The imaginary part of the relative permittivity can be ob- 
tained as the product of the relative dielectric constant 
and the loss tangent. 

Lundien found the behavior of the dielectric constant 
to be very similar for all soils tested. The relative di- 
electric constant increases with increasing water content 
from a minimum of approximately 3. Over an intermediate 
range of moisture levels, the relative dielectric constant 
of sand is generally about 151 greater than that of clay. 

Of course there are numerous cases which fall within the 
extremes , 

Lundien' s measurements show no relation of the dielec- 
tric constant to the compaction of the samples. Various 
samples were prepared with three different compactive 
forces; S.74, 11.83, and 18.47 N/cm^. The lack of varia- 
tion of the data due to differing compaction seems to in- 
dicate that the dielectric constant is more dependent on 
the volumetric water content [grams per cubic centimeter) 
than on the percent moisture by weight. 

Considering the data obtained for all samples, Lundien 
formulated an average relation between water content and 
relative dielectric constant as follows: 

VWC = ~ 
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3 

Inhere VWC = volumetric water content, g/cm 

€ ’ - relative dielectric constant 
r 

The range of frequencies over which Lundien measured 
the dielectric constant was too small to provide conclusive 
information about the effects of frequency on the dielec- 
tric constant. However, comparison of Lundien' s results 
with data from a previous study [6] indicates that the 
dielectric constant of sand and silt is approximately con- 
stant over the range of frequencies from 0.3 GHz to 1.5 
GHz. The permittivity of various soils has been measured 
at higher frequencies by other investigators, as discussed' 
below. 

Measurements by Wiebe 

The relative permittivity of nine types of Texas soil 
was measured by Wiebe [7], Several samples of each soil 
were prepared with different moisture contents and compac- 
tion. The results were very similar to those obtained by 
Lundien. 

Wiebe investigated two techniques for measuring the 
relative permittivity. Figure II-2 illustrates the ar- 
rangement of equipment for the "free space” technique in 

which the sample to be measured is placed between two an- 

\ 

tennas. For the other technique, . the iantennas : are. removed 
and a section of waveguide inserted. A portion of the 




Figure II- 2. Wiebe's arrangement for measuring 
relative permittivity of soil [7]. 
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waveguide is removable and is used to contain the sample. 

The procedure is essentially the same for both tech- 
niques. With no sample in place, the phase and attenuation 
are adjusted to give a convenient reference level on the 
indicator. When the sample is inserted, the attenuator 
and phase shifter are readjusted to give the same reading 
on the indicator. The difference in attenuation and phase 
from the initial settings are the attenuation and phase 
shift due to the sample. By varying the length of the 
sample, the attenuation and phase shift can be plotted as 
functions of sample length. The slopes of these curves 
give the attenuation constant a and the phase constant 3 . 
The relative permittivity of the sample is determined from 
one of the following sets of equations. For the free space 
technique, 



(II-23a) 


Cll-23b) 


where 3o is the phase constant of free space, and for the 
waveguide method. 


€. = 



t 


(II-24a) 
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where K is the wave number evaluated at the cut-off wave- 
c 

length of the guide. 

Wiebe made measurements with both techniques at 9.0 
GHz to evaluate the quality of the measurements. Satisfied 
that both methods gave good results, Wiebe chose to make 
further measurements with the waveguide method at a fre- 
quency of 10.625 GHz. 

Wiebe' s measurements are generally consistent with 
those of Lundien. The behavior of the permittivity as a 
function of moisture content is essentially the same for 
all soils tested. In addition, the permittivity of sand 
is found to be 15% to 20% greater than that of clay over an 
intermediate range of moisture. 

By measuring the samples immediately after wetting and 
after 24 hours of curing, Wiebe was able to show the change 
of permittivity due to adsorption of water. The permittiv- 
ity is highest immediately after wetting, when most of the 
water is still in a free state. The effective permittivity 
of bound water is less than that of free water. Therefore, 
as more water molecules adhere to the soil, the average 
permittivity of the water decreases, thereby decreasing the 
permittivity of the soil-water mixture. 

Some of Wiebe 's data show more dependence on 
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compaction than Lundien's data. The apparent difference 
may be due to the fact that Wiebe used percent moisture 
content (by weight) as a measure of water content of a sam- 
ple rather than volumetric water content, as used by 
Lundien. 

Ground-Based Measurement of the Apparent 
Temperature of Smooth Sand 

Having knowledge of the permittivity of sand, Richer- 
son [8] performed an experiment to test the accuracy of 
Peake's model for apparent temperature. A specialized form 
of Peake's model was used to predict the apparent tempera- 
ture of a smooth bed of sand for various moisture contents. 
Data were obtained with a 31.4 GHz radiometer mounted on a 
small tower adjacent to the test surface. 

Richerson modified Peake's model by inserting the re- 
lation 

>iiH = 

for the differential scattering coefficient of a smooth 
surface, into (II-3). It follows that the albedo of the 
smooth surface is 

“ [Ri(^o)| (II -25) 

where R^[0o) is the Fresnel reflection coefficient. In 
addition, Richerson added a third term to (II-ll) to 
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account for radiation of the atmosphere between the surface 
and the radiometer. 

The measured apparent temperatures were not in good 
agreement with predictions. The moisture content of sur- 
faces for which results were reported varied from 7.56% to 
13.2%. Over this range of moisture there was no signifi- 
cant change of apparent temperature. 

Several reasons for the discrepency between theoreti- 
cal and measured apparent temperatures can be cited. By 
considering the effects of a non-ideal antenna pattern, 
Richerson was able to predict smaller variations of appar- 
ent temperature due to moisture. Furthermore, Richerson 
pointed out that the short wavelength (9.55 mm) sensor is 
very sensitive to surface characteristics. Not only did 
roughness effect the measurements, but drying of the soil 
at the surface caused significant changes in the apparent 
temperature. 

Experiments similar to Richerson’ s have been performed 
by other investigators with active systems. The results 
of the experiment described below show more clearly the 
effect of surface roughness at various frequencies. 

Ground-Based Measurement of the Reflectivity 

of Moist Soil 


Stiles, et al . [9], measured the power reflectivity of 
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some test surfaces with various moisture contents to deter- 
mine whether moisture content could be monitored by radar. 
Measurements were made with a swept-frequency system with a 
frequency range of 4.0 to 26.5 GHz., The results show a 
definite relation between moisture content and reflectiv- 
ity. However, measurements of rough surfaces indicate an 
"apparent moisture" less than the actual moisture content. 

Data obtained from smooth sand at 6.0 GHz show a def- 
inite increase of reflectivity with increasing moisture 
content. The relation is apparently linear. A similar 
relation is evident at 26 GHz, but the reflection coeffi- 
cient does not increase as rapidly (as a function of mois- 
ture) as at the lower frequency. 

The reflectivity of a rough surface is lower than that 
of a smooth surface for all frequencies. However, the dif- 
ference between reflectivity of smooth and rough surfaces 
is much greater at higher frequencies. The roughest sur- 
face tested by Stiles had a reflection coefficient of 0.28 
at 8 GHz and 0.04 at 26 GHz, while the reflection coeffi- 
cient of the smooth surface was almost constant at 0.36 
over the same frequency range. These data emphasize the 
fact that there is no absolute scale for roughness. Any 
rough surface may be considered smooth for sufficiently 
long wavelengths. 

The results of the, experiment .prove the possibility of 
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monitoring the moisture content of smooth surfaces with 
radar. However, the difficulty of remotely measuring the 
moisture content of an arbitrarily rough surface is also 
apparent. 

Airborne Measurements of Apparent Temperature 

To better determine the feasibility of measuring soil 
moisture content with passive microwave sensors, personnel 
of NASA/Goddard Space Flight Center have conducted flights 
over selected agricultural test sites to gather microwave 
radiometer data. Ground support data for each flight in- 
cluded the moisture content of several soil samples from 
each field along the flight line. The data from two such 
flights have been analyzed at Texas A§M University to 
determine the correlation of apparent temperature with 
moisture content. 

Apparent Temperatures of Bare Fields 

Data obtained from Weslaco, Texas, were analyzed by 
Jean [10]. The data included apparent temperatures mea- 
sured at 1.42, 2.69, 4.99, and 10.69 GHz. The fields at 
Weslaco were rough but not vegetated. The moisture content 
varied from 6.71 to 35%. 

To quantify the correlation of apparent temperature 
with moisture content, Jean computed a rank difference 



29 


coefficient of correlation for each set of data. This 
coefficient is defined by 

P -l) 

where = rank of the element in data set {X} 

= rank of the element in data set {Y} 

n = number of pairs of data {x^, y^} 

Considering the values of the correlation coefficient for 
each data set, the 2.69 GHz data showed the best correla- 
tion and the 1,42 GHz data the next best. However, none 
of these showed a strong correlation. 

To illustrate graphically the correlation of apparent 
temperature with moisture, Jean plotted the average appar- 
ent temperature of each field as a function of the average 
moisture content. Figure II-3 is an example of the re- 
sults. The data are from the horizontal channel of the 
1.42 GHz sensor. The straight line was fit by a least 
squares technique. The slope of the lines is -1.2°K per 
percent moisture. 

The data presented by Jean clearly indicate the poten- 
tial for remote measurement of soil moisture in bare 
fields. However, similar results have not yet been obtain- 
ed for vegetated fields. 
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Apparent Temperatures of Vegetated Fields 

An analysis similar to Jean’s has been performed by 
Kroll [11] on data obtained from Chickasha, Oklahoma. The 
data include apparent temperatures measured at 1.42, 4.99, 
10.69, and 37 GHz. 

The surface conditions of the test fields at Chickasha 
were quite different from those at Weslaco. The fields at 
Chickasha were relatively flat. Furthermore, most were 
covered with young wheat varying in height up to about 25 
cm. 

Figure II- 4 is a plot of apparent temperature versus 
moisture content for the 1.42 GHz data. There is no clear 
relation between moisture and apparent temperature. This 
may be due partly to the limited range of moisture avail- 
able. However, it may be that the effects of the vegeta- 
tion suppress the effects of varying moisture in the soil. 

The theories and experimental results presented in 
this chapter are not adequate to describe emission and 
scattering from soil beneath a plant canopy. The diffi- 
culty of predicting quantitatively the apparent temperature 
of even a smooth surface of known permittivity has been 
shown . However, other experimental data indicate that the 
apparent temperature of bare soil is indeed a function of 
moisture content. An attempt to predict the apparent tem- 
perature of some simple vegetated surfaces is discussed in 



Percent Soil Moisture 


Figure II-4. Average apparent temperature at 1.42 
GHz for fields at Chickasha, Oklahoma 
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the following chapter . 
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CHAPTER III 

APPARENT TEMPERATURE OF VEGETATED TERRAIN 


The apparent temperature of any object is a measure of 
the thermal radiation emanating from the body in a parti- 
cular direction with a particular polarization. For micro- 
wave frequencies the apparent temperature is related to the 
radiant power through the Ray leigh- Jeans approximation of 
Planck's Blackbody Radiation Law [2], 


T — ^ T wgf 


(in-i) 


where k * Boltzmann's constant 

T = temperature in degrees Kelvin 
Af = a narrow band of frequencies 
X = a wavelength of center frequency 

Only for a blackbody, that is, an ideal absorber, does 
(III-l) give the relation between the thermal radiation and 
the actual temperature of the body. For any real object 
the temperature T of (III-l) is the apparent temperature, 
but is not the actual temperature of the body. 

The apparent temperature of vegetated terrain results 
from contributions of radiation from several sources. The 
most significant contribution is the thermal emission from 
the ground. In general this contribution includes 



35 


radiation from both the soil and the vegetation. A second 
contribution is from part of the diffuse radiation of the 
atmosphere which is scattered by the ground toward the sen- 
sor. Scattering may be due to both soil and plants. An- 
other contribution may be due to the scattering of radia- 
tion from quasi-point sources into the antenna beam. This 
contribution is usually negligible [2]. The final contri- 
bution to the apparent temperature of the ground is from 
the radiation of the atmosphere between the ground and the 
sensor. 

The purpose of this chapter is to determine the vari- 
ation of the apparent temperature of vegetated terrain as 
a function of soil moisture content. Calculation of the 
apparent temperature is facilitated by considering only 
thermal emission of the soil and vegetation. Atmospheric 
effects are discussed in Appendix B. 

Elements of the Transmission Problem 

The determination of the apparent temperature is pri- 
marily a problem in propagation through dielectric layers.. 
Part of the thermal radiation which originates within the . 
soil is transmitted across the surface of the soil into the 
canopy. The transmitted radiation then propagates through 
the canopy and into the atmosphere. 

The most useful descriptors of the dielectric media 
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are their permittivities. These are discussed below. 
Permittivity of the Soil 

As noted in Chapter II, the permittivities of various 
types of soil have been measured at several microwave fre- 
quencies [5], [7]. For most soils a curve of relative per- 
mittivity as a function of moisture content is very similar 
to Figure III-l. Although the permittivity is a function 
of frequency, data presented by Wiebe [7] indicate that the 
change in permittivity over the frequency range from 1.42 
GHz to 31.4 GHz is no more than might be caused by a varia- 
tion in soil type. Therefore, the curves of Figure III-l 
may be assumed to represent the relative permittivity of 
some particular soil at any particular frequency in that 
range. 

Permittivity of Vegetation 

It is not surprising that the permittivity of vegeta- 
tion is also a function of water content, since plants are 
composed largely of water. From observation of X-band 
data, Peake and Oliver [12] have formulated the following 
approximation for the permittivity of vegetation: 

+ j (in-2) 
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where F = fraction of water by weight in the plant 
= permittivity of water 

This approximation is valid only for frequencies above 
1 GHz, in which case the ionic conductivity of dissolved 
salts is negligible [12]. 

Water has been found by experiment to exhibit Debye 
type relaxation [13], Debye relaxation is defined by the 
exponential behavior of the displacement current in a di- 
electric to which an electric field is suddenly applied 
[14]. This behavior is expressed as 

D(t >0) = + (e, - ej(|- e ^)E 

where = instantaneous dielectric constant 
Eg = Static dielectric constant 
T = relaxation time 


If an electric field having sinusoidal time variation is 
applied to such a dielectric, the complex permittivity is ' 


€ 


- e 


oo 


*oo 


(III-4) 


I +j27rfT 

For application to water at microwave frequencies it 
-Onvcnieiit to convert ZttT to an equivalent frequency. 


Thus, (III-4) may be written 


= € 






oo 


1 + 


m.) 


(III-S) 



39 


The high frequency dielectric constant is constant, 
but the static dielectric constant and the relaxation 
frequency f^ are functions of temperature. There is some 
disagreement about the value of but 5,5 is acceptable 
[13], The static dielectric constant is 

= 87.7 - 0.4(T- 273) 

where T is the temperature in degrees Kelvin. For tempera- 
tures between 273°K and 303°K, is approximated by 

f^ =• 9.0 + 0.40S(T-273) GHz 

Using these parameters in (III-5) yields the value of 
for (III-2). 

Permittivity of a Dielectric Mixture 

The permittivity given by (III-2) is that of the plant 
material only. The canopy above the soil is a mixture of 
vegetation and air. Since the canopy is treated as a sin- 
gle medium in the analysis which follows, it is necessary 
to determine the permittivity of the plant-air mixture. 

A model applicable to the present case is the Weiner 
model for a dielectric mixture as presented by Evans [15]. 
The relative permittivity of ;the' mixture : is .related. to the 
relative permittivities of the constituent media through 
the formula. 
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^2 ■/ 
€2 + u 


CIII-6) 


in which p is the fraction of the total volume occupied by 
medium 1. 

The parameter u in (III-6) is called the Formzahl. 

This parameter describes the dispersion of one medium with- 
in the other. The value of u may be any non-negative num- 
ber. The significance of several values of the Formzahl is 
illustrated in Figure III-2. It seems very difficult to 
select an appropriate value of u for a plant canopy. How- 
ever. Figure III-3 illustrates that the selection of u is 
not critical in the range u = 10 to u = 25, since the rela- 
tive dielectric constant varies by only about 15% around 
the mean value. 

The derivation of an explicit expression for c is 
facilitated by assuming that medium 2 is air with relative 
permittivity equal to unity. Thus, the last term of 
(III-6) vanishes, leaving 

IjhzI - p iiJii 
e^ + u r e,+u 


which leads to 
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Direction of Electric Field 



Figure III-2. Various types of mixture represented 
by the Formzahl. (After Evans (7]) 
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e,(l + pu) + u(l-p) 
e,(|-p) + p + U 


(III-7) 


where is the relative permittivity of the vegetation, ^ 
and p is the fraction of the canopy volume which is occu- 
pied by vegetation. 

Smooth, Uniformly Vegetated Surface 


The simplest model for vegetated ground is a smooth 
surface covered by a continuous canopy of uniform height. 
The soil is assumed homogeneous with permittivity as illus- 
trated in Figure III-l (p. 37). The plant canopy is repre- 
sented by a homogeneous dielectric layer with permittivity 
as given by (111-7). The atmosphere is assumed lossless 
and nonradiating, 

Threer.processes determine the apparent temperature of 
the vegetated ground. First a portion of the thermal radi- 
ation incident on the surface from below is transmitted 
across the boundary toward the sensor. The transmitted 
field experiences attenuation as it propagates through the 
canopy. Finally, the thermal radiation from the canopy 
augments the radiation from the surface which propagates 
through the canopy. 
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Transmission Coefficient 

The power radiated into the plant canopy is equal to 
the power incident on the surface from below times the 
transmission coefficient of the surface. The transmission 
coefficient is defined by 

T^ = Power transmitted with polarization i 
Power incident with polarization i 

and can be expressed in terms of the permittivities of the 
two media. For this derivation the second medium is as- 
sumed to be air, since it is usually air which forms the 
boundary with the soil. 

The relation between the incident and transmitted 
fields is shown in Figure III-4 for vertical and horizontal 
polarizations. The angles 0 and ijj are related by Snell's 
law of refraction [16], 

Vj -- (III-8) 

where v^ = velocity of propagation in air 
v^ = velocity of propagation in soil 

The velocity of propagation in air may be assumed equal to 
the free space velocity of light. The velocity of propa- 
gation in soil is 
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Figure III- 4. Orientation of incident and trans- 
mitted fields. 
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where e' = real part of permittivity 

e" = imaginary part of permittivity 

The incident and transmitted power densities (watts/ 
square meter) are determined by Poynting's theorem [17], 

For the case of horizontal polarization [Figure III-4(a),, 
p. 4S] the incident fields are 

E, = E, X 

H, = (f - (f 

where n-|^ = is the intrinsic impedance of the soil. 

Performing the cross multiplication leads to 

T + zc^V')] 

CIII-9) 

for the incident power density. The power incident upon a 
surface element dS is 
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O^S (III-K 

The transmitted fields are 

E,= r,E.x 

H, = ^(hjCo^O - ZxUmO) 

where is the transmission coefficient for the horizon-, 
tally polarized electric field. 

_ 2 \fiT 

CO^O Cin-lj 


Applying Poynting's theorem yields 
or 



for the transnvitted power density. The power transmitted 
from a surface element dS is 

Pavj " x|E,n7;|^(^)c(So2^ (in-13) 

The desired transmission coefficient is simply the 
ratio of (III-13) to (III-IO). 
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A similar analysis for the case of vertical polarization 
[Figure III-4(b), p. 45] leads to 



\ iJ 


where 


2 

Vi7 + Vi7 oo<i9 


(III-15) 


(III-16) 


Attenuation by Canopy 

Since the canopy is represented as a lossy dielectric, 
an attenuation constant can be determined for the canopy. 
The attenuation constant a is defined as the real part of 
the complex propagation constant of the electric field. 

r = a + j/3 

Noticing that the propagation constant is also written [17] 
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a solution for the attenuation constant in terms of the 
permittivity can be obtained. The result is 



(III-17) 


Since the ratio a/to is equal to the imaginary part of the 
permittivity, (III-17) may be written in a more convenient 


form. 


a - 



where = real part of permittivity 

e” = imaginary part of permittivity 


From (III-9) it is clear that the average power pro- 
pagated is proportional to the square of the magnitude of 
the electric field. The magnitude of a field having pro- 
pagated a distance D through a medium with attenuation 
constant a is 


E 



where is the initial field strength. therefore, the 
power at the same point is 


P 


D -2CiD 

R. e 


The power radiated from the soil may be. represented by 
its apparent temperature. The apparent temperature just 
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above the surface is equal to the product of the ground 
temperature, the emissivity of the soil, and the transmis- 
sion coefficient. Measurements of the apparent temperature 
of smooth sand [7] indicate that the emissivity is greater 
than 0.95. However, for a qualitative analysis such as the 
present, any value of emissivity may be. assumed. Consider- 
ing the attenuation by the canopy, the contribution of the 
soil to the total apparent temperature is 



_ T -2o.WAtcB 


(III-18) 


where T = ground temperature 

o 

= emissivity of soil 

= transmission coefficient for polarization i 
H * canopy height 


Radiation from Canopy 

The canopy not only decreases the apparent temperature 
of the soil, but also contributes to the total apparent 
temperature through thermal emission. Although most of the 
radiation is emitted by the leaves of thecplants, it is 
assumed that all points within the canopy emit radiation 
equally . 

The contribution of the canopy to the total apparent 
temperature may be derived from a general expression for 
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the apparent temperature of a dielectric layer of thickness 
H. The expression is [18] 

H o(z] 

\ = f McOdU ciii-19) 

0 

where T(h) = temperature of medium at height h 
£^(h] = differential emission coefficient 
a(z) = power attenuation constant at height z 


Assuming that the temperature, attenuation constant, and 
emittance are constant within the canopy, (III-19) becomes 


( k^opy ^ 

where a is defined by (III-17). 
leads to 


~2aAec&(H -li) 

S dh 

Performing the integration 



anopy 
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(IIIr20) 


According to Kirchhoff*s law [2], if the radiating 
layer is in thermal equilibrium with its environment , . the 
differential emission coefficient (radiation per unit 
length) is equal to the absorption coefficient, expressed 
in (III- 20) as 2a. As stated in Chapter II, the assumption 
of thermal equilibrium is often violated by a plant canopy. 
It may be stated quite generally, however, that the emis- 
sion coefficient is equal to the product of 2a and an 
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energy transfer factor f. Thus, (III-20) may be written 





-2aH>dec<9 j 


CIII-21) 


where f < 1 represents a gain of energy by the canopy, and 
f > 1 represents loss of energy by the canopy. 


Theoretical Predictions 


For reference, the apparent temperature of a smooth 
soil surface is illustrated in Figure III-5. The apparent 
temperatures for horizontal and vertical polarizations are 
plotted as functions of refraction angle for two values of 
soil moisture content. With vertical polarization the ap- 
parent temperature has a characteristic peak at the Brew- 
ster angle. It is clear that for angles less than the 
Brewster angle, an increase in moisture causes a decrease 
in the apparent temperature. 

The total apparent temperature of the vegetated sur- 


face is obtained by combining (III-18) and (III-21). 





-2aHxlt<:e 
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(III-22) 


) 


Apparent temperatures computed from (III-22) are 
plotted as functions of soil moisture content in Figure 
III-6. The vegetation was assumed to occupy 5% of the 
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total volume of the canopy. The height of the canopy was 
varied from 25 cm to 100 cm. It is clear that increasing 
the height of the canopy decreases the sensitivity of the 
apparent temperature to variations in soil moisture. 

Increasing the density of vegetation has the same 
effect as increasing height, as illustrated in Figure 
III- 7. To compute the data of this figure, the canopy 
height was assumed to be 50 cm, and the percentage of can- 
opy volume occupied by vegetation was varied from II to 
101. From observation of fields of grain sorghum and cot- 
ton, it has been determined that the density of these crops 
may be less than II. However, depending on the method of 
planting and amount of growth, the vegetation density may. 
be 31 to 51. It can be seen in Figure III-7 that as the 
density of vegetation increases, the apparent temperature 
increases and becomes less sensitive to variations of soil 
mois ture , 

The significance of the contribution of the canopy to 
the apparent temperature also depends on the frequency of 
radiation considered, as illustrated in Figure III- 8. A 
canopy height of 50 cm is assumed. At 3 GHz the sensitiv- 
ity of apparent temperature to soil moisture is approxi- 
mately 1.5°C/percent moisture. At 5 GHz the sensitivity is 
about 0 . 5°C/percent moisture, and at 7 GHz the dependence 
on moisture is almost negligible. 



Apparent Temperature (°K) 
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Row Crops 

An expression for the apparent temperature of row 
crops can be developed using the parameters derived for the 
case of uniform vegetation. More variables must be in- 
volved, since the apparent temperature of row crops depends 
not only on the height and density of vegetation, but also 
on the distance between rows and the angle from which the 
ground is observed. 

The geometry of the problem is illustrated in Figure 
III-9. For simplicity the rows are assumed to be rectan- 
gular in cross section. The medium which comprises each 
row is a homogeneous dielectric with permittivity defined 
by [III-2) and (III-7), The surface of the soil is assumed 
to be an infinite plane. The direction from which the 
terrain is observed is defined by a refraction angle, 6, 
and an azimuth angle, , both of which assume only values 
less than or equal to it/ 2. The direction 9 = 0 is normal 
to the surface, and = 0 is perpendicular to the rows. 

Emission Perpendicular to Rows 

Before the general form of the apparent temperature 
of row crops is derived, the special case of tj) = 0 is con- 
sidered. 

A cross-sectional view of a few rows is given in Fig- 
ure III- 10. The height and width of each row are denoted 
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H and W, respectively. The width of the non- vegetated 
space between rows is denoted A. A dimension B = H tan 9 
is also defined. 

The length CA + W) may be considered a unit length or 
a period. The power radiated into the direction (0, by 
any region of length (A + W) is equal to that radiated into 
the same direction by any other region of equal length. 
Therefore, the apparent temperature of a large area is 
equal to the average apparent temperature of the region 
specified in Figure Ill-lO(a). 

If B < A as in Figure Ill-lO(a), only a portion of the 
radiated power experiences attenuation. The distance tra- 
versed through the canopy depends upon the point on the 
surface from which the radiation emerges. The average 
length of the path through the canopy is equal to the 
cross-sectional area of the canopy divided by the projec- 
tion of (B + W) into a plane normal to the refraction di- 
rection. 

HW 

= (B 


where d is the average distance traversed through the can- 
opy in the direction (0, ” 0)* 

The total apparent temperature is the average of the 
contributions from the bare soil and the vegetated area. 
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In general, 


(A-B)lU„.i-(B-fW)(U.. 
A + W 


(III-24) 


The contribution of the bare soil is simply 

where T = temperature of ground 
& 

e - soil emissivity 

= transmission coefficient for polarization i 

The contribution of the region in which the soil is covered 
or shadowed by vegetation is given by ( 111 - 223 , which is 
repeated here. 






Ves 


(III-26) 


+ TJ(/ 


in which the attenuation distance Hsec 6 has been replaced 
by d. 

A somewhat different expression is obtained for the 
case in which B > A, In such a case, all of the soil is 
either covered or shadowed by vegetation. Consequently, 
the total apparent temperature is given by (III-26), and 
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the average attenuation distance is 


HW 

(a VJ)co^9 


(III-27) 


Emission at Any Azimuth Angle 


The formulas which express the apparent temperature 
of row crops for the special case of = 0 can easily be 
modified for the general case, 0 £ (J> j< n/2. 

Figure III-ll illustrates the geometry for propagation 
through the canopy at an arbitrary angle. The dimension B 
is defined generally as 

B = (III-28) 


The projection of the canopy cross section onto the plane 
of observation is HWsec (p. As before, the average attenua- 
tion path is expressed as this cross-sectional area divided 
by the projection of (B + W) onto a plane normal to the 
direction of observation. 


H 

(B + W)xt2C.0 ocxt$ 
HW 

(B +\f\/)oo<0 


(III-29) 
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which is identical to (III-23). 

It follows that the other equations developed for the 
special case <J> = 0 also are true for the general case. 
Those expressions which apply when A > B and those which 
apply when A < B may be combined into a single expression, 
by defining a special function as follows. 


L(A,B) = B when B < A 
L(A,B) = A when A ^ B 

The apparent temperature of row crops may now be expressed 
as 

^ _ [a-L(A,B)](TJ,„/[UA,B)-w1(tJ„. 

T,i ^ (H.-SO) 


where the contribution of the bare soil is given by (III- 
25), the contribution of the vegetated soil is given by 
(III- 26), and the average attenuation distance is 


HW 

[L(A,B) +• V^]oo^9 


(III-31) 


Theoretical Predictions 

The dimensions selected to describe the theoretical 
row crops approximate the dimensions of certain actual 
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crops. Fields o£ young cotton and grain sorghum have been 
observed to determine the height-width ratio of these 
plants. Since the standard row spacing of these crops is 
40 inches, the dimension (A + W) has been set at one meter. 

Apparent temperatures computed from (III-30) are 
plotted in Figures III-12(a) through [c) as functions of 
soil moisture. In general, the apparent temperature de- 
creases as the moisture content increases. The data for 
both polarizations exhibit the same behavior, although the 
vertical polarization data are consistently higher than the 
horizontal data. However, as the amount of plant cover 
increases, the sensitivity to changes of soil moisture 
decreases, as well as the difference between horizontal and 
vertical polarization measurements. 

Figures III“13(a) and (b) illustrate the dependence of 
the apparent temperature on the angle from which the ter- 
rain is observed. According to the model developed, vege- 
tation causes a net increase in the apparent temperature of 
the ground. Therefore, the apparent temperature of row 
crops is greatest when measured from large refraction 
angles and small azimuth angles, in which cases the bare 
soil between rows is shadowed by the plants. The apparent 
temperature decreases as the azimuth angle is increased, 
reaching a minimum at (J) = 90. 
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Figure III-12 (b) . Apparent temperature of row crops 

as function of soil moisture. 
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A = W = 50 cm 
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Figure III-13 (a). Apparent temperature of row crops 

as function of azimuth angle. 
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Figure III-13 (b) . Apparent temperature of row crops 

as function of azimuth angle. 
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Emission from a Rough Surface 

The models developed for the apparent temperature of 
vegetated terrain are ideal, since smooth surfaces have 
been assumed. No such surface is encountered in practice. 
However, for cases in which surface variations are small 
with respect to a wavelength of the radiation being mea- 
sured, the smooth surface model may be a reasonable approx- 
imation. It is desirable to have a model applicable to a 
rough surface. 

A model for emission of * thermal . radiation . from a rough 
surface has been developed by Johnson [19]. Johnson's 
model expresses the radiation refracted into a direction 
(0 , 4) ) in terms of the radiation incident from (0., d>.) 
and a "coefficient of effective area". The results of 
Johnson's derivation are presented below and modified for 
application to the measurement of apparent temperature. 

Transmissive Scattering Coefficients 

The emission of radiation from a rough surface is 
referred to by Johnson as "transmissive scattering" since 
the process is similar to ordinary scattering. Electro- 
magnetic energy incident on the surface from a direction 
(®i> is refracted into various directions depending 
on the characteristics of the surface. The intensities of 
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the fields refracted into a particular direction (6^, 
are given by the transmissive scattering coefficients, 


defined as follows: 




Ue. >■) 

E«l^i Al 



HhI^. a) 

H,(^, M 

mA> 

H,(^, ,'^i) 

HMAI 


The subscripts H and V denote horizontal and vertical 
polarization with respect to the mean .coordinate system. 

For the transmission coefficients the first subscript 
indicates either electric or magnetic fields, the second 
indicates polarization of the scattered field, and the 
third indicates the polarization of the incident field. 

The development which follows is based on the assump-. 
tion that each element of surface area has an "effective 
area" oriented in such a way as to cause specular refrac- 
tion from 0^) into (0^, (I>g). 
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Localized Parameters 

Each element of surface area is described by two nor- 
mals, one internal and one external. Each of these normals 
is defined by a zenith angle and an azimuth angle with re- 
spect to the mean coordinate system. Internal and external 
coordinate systems may be selected so that both normals are 
defined by a single directional pair (0^, ‘ 

The direction of the effective normal is implied by 
the incident and scattering directions. The normal lies 
in the plane formed by the incident and scattering direc- 
tions. Furthermore, the normal must be such that the local 
angle of incidence, and the local angle of refraction, 

6gn» related by Snell*s law of refraction, 

where is the velocity of propagation in the lower 

medium, and is the velocity of propagation in the upper 
medium. 

The local angles of incidence and refraction are re- 
lated to other parameters by the Law of Cosines for sides 
of oblique spherical triangles, 

(III -32a) 
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An incident field polarized either vertically or hori- 
zontally with respect to the mean coordinates appears to 
have both vertical and horizontal components with respect 
to the local coordinates. The locally polarized components 
can be obtained by rotating the incident polarization coor- 
dinates about the line of incidence. The angle of rotation 
is determined from the Law of Cosines for sides of oblique 
spherical triangles, 

ocri - CokiB C<ki9; 

Co^a, - — - — — (in-33) 

The mean polarization of the scattered fields is obtainable 
by a similar rotation of local polarization coordinates 
about the line of refraction. The angle of rotation is 
determined by 

ca!^ 9„ -OKi9s^c^9s 

C/0<tay, = (III-34) 


Effective Surface Area 


Each element of surface area has a corresponding 
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projection on the mean coordinate plane. Except for a 
perfectly smooth surface, the actual surface area is great- 
er than the projected mean area. The "surface roughness 
factor" is defined as the actual surface area per unit 
projected area, 


8 


dA 


• 8 i I 


(III-35) 


The refraction of radiation from (0^, 4)^) into 
(0 , ) implies an effective area oriented so as to cause 

w w 

specular refraction. A "coefficient of effective area" is 
defined to relate the effective area to the actual surface 
area. 

d A^ff = /3dS ; 0 i P < \ cin-36) 


where 3 is the coefficient of effective area. 

The roughness of the surface may be described by the 
distribution of the normals. The zenith and azimuth angles 
are treated as random variables with joint probability 
density 4’j^) • Assuming that the azimuth angle is 

uniformly distributed and independent of the zenith angle, 
the joint density may be expressed as 

where is the marginal density of the zenith angle 
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of the normal. 

The coefficient of effective area and the surface 
roughness factor are related to the statistics of the sur- 
face normals. The probability that the normal is oriented 
in the proper direction, within an interval dS , d(J) is 
(l/27r)fQCe )d0 d(J). Since a solid angle bounded by d0 , d4> 
is defined as 

otD. - dO d4^ 


the probability that the surface normal is oriented within 

dQ is (l/2TTsine) (0^) dj}. The coefficient of effective 

b n 

area is 





(in-38) 


which is the fraction of the total surface oriented within 
a unit solid angle about a particular direction. The pro- 
jection of this area on the mean plane is 3cos0^. The 
mean projected area per unit surface area is obtained by 
integrating over the hemisphere. 

da 

dS J J 27 t 


(III-39) 
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From (III-35) and (III-39) the surface roughness v factor is 



(III-40) 


From (III-35) and (III-36) the effective area may be 
expressed generally as 


dA^ff = /GS dA 


Substituting (III-38) and (III-40), 



CIII-41) 


which expresses the effective area causing specular re- 
fraction from (e^, <f)^) into (0^, as a function of 
surface parameters. 


Apparent Temperature of the Surface 

As stated earlier in this chapter, the apparent tem- 
perature of any object is a measure of the radiant power 
emanating from the body. Johnson expresses the polarized 
spectral radiance of a rough surface in terms of the in- 
ternal radiation, the transmission coefficients, and the 
surface characteristics. Correcting for an apparent trans- 
position error, Johnson's expression is [19, Eq. IV-43] 
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■ CC<Ld-. 


7^ 7“ I 

EVV HVV ■*’V 

^EVH ^HVH 


fill 


(ni-42) 


N„ = 


S / j 


^EHV ^HHV^v) 


BHV ' HHV 

where Ny, = polarized component of radiated power 

ly, Ipj = polarized component of internal incident 
power 


The electric field transmission coefficients are 




^Evv " C(rta^ + 7^^^ /ci4>naj 


^£VH " ^<Uoria^ 


^EHV V '*’ CLj C/O^CL^ 

where and Tcu are the Fresnel transmission coeffi- 
b V nn 

cients for vertical and horizontal polarization. The mag- 
netic field transmission coefficients are given by the same 
expressions, changing each subscript E to H. 
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The radiance expressed by (III-42) is the power re- 
fracted into the direction (0g» due to power incident 
from a single direction (0^, . The various orientations 

of surface elements allow radiation from other incident 
directions to be refracted into (©«.» 4’c)* Assuming that 
the vertical and horizontal components of the incident 
radiation are equivalent, the apparent temperature of a 
rough surface without vegetation is 




(in-43) 

Cxw V Eiv 'hiv 


where the subscript I denotes either vertical or horizontal 
polari zation . 

The total apparent temperature of a vegetated rough 
surface is obtained from the equations derived earlier in 
the chapter by inserting (III-43) for the contribution of 
the soil. 



CHAPTER IV 


SCATTERING COEFFICIENT OF VEGETATED TERRAIN 

Whenever electromagnetic energy is incident upon some 
object, a portion of that energy is reflected, or reradi- 
ated from the body. If the surface is very smooth the re- 
flected energy due to an incident plane wave may be concen- 
trated into a single direction. Such reflection is termed 
specular . The energy reflected from a rough surface is 
generally scattered into various directions. Such scat- 
tering is termed diffuse . In general, the reradiated 
energy from any object includes both specular and diffuse 
components.. 

The radar cross section of an object is defined as the 

area required to intercept an amount of power which, when 

scattered isotropically, produces the same power density at 

the receiver as the target [20]. The power density (watts/ 
2 

m ) of the incident radiation is 



where is the impedance of the atmosphere. The cross- 
sectional area of an ideal, lossless, isotropically scat- 
tering target is denoted a. Therefore, the power scattered 
from the ideal target is 
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and the power density at the radar receiver, a distance R 
from the target is 

rp = o-IEil" = JlJ" 

where is the electric field intensity at a distance R 
from the target. Thus, the radar cross section of an arbi- 
trary target is defined to be 

O' = AttR^ 

regardless of the actual cross-sectional area. 

An extended target which is too large to fall within 
the radar antenna beam is best described by the radar cross 
section per unit area, also known as the scattering coeffi- 
cient. This parameter is usually denoted a®. It may be 
defined as the radar cross section of the target divided 
by the area of illuminated surface. 

In general, scattering from vegetated terrain includes 
scattering from both the soil and the plant canopy. How- 
ever, one of the two components may have a predominant 
effect on the total scattering coefficient. In this chap- 
ter models are developed for the backs cattering coeffi- 
cients of certain types of vegetated terrain to determine 
the dependence of the backs catter coefficient on the mois- 
ture content of the soil. 


(IV- 1) 
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Elements of the Scattering Problem 

Energy incident upon the vegetated surface is gener- 
ally scattered by both the canopy and the soil. Further- 
more, that portion of the incident field which eventually 
reaches the soil, experiences additional attenuation due 
to absorption within the canopy. Similarly, the radiation 
scattered by the soil surface is attenuated as it propa- 
gates upward through the canopy. 

The soil and the canopy are described by their per- 
mittivities as in Chapter III, The soil is assumed to 
have the relative permittivity illustrated in Figure III-l 
(p. 37) as a function of moisture content. The permittiv- 
ity of vegetation, computed from (III-2) is used to deter- 
mine the equivalent permittivity of the canopy according to 
(in-7). 

Smooth, Uniformly Vegetated Surface 

The first type of terrain considered is a smooth sur- 
face covered by a uniform layer of vegetation. For the 
purpose of computing the attenuation due to the vegetation, 
the canopy is assumed to be a homogeneous dielectric layer. 
The backscatter coefficient of the soil is derived by the 
physical optics method. Scattering from the canopy is 
accounted for by a modified form of Peake ’ s :model for scat- 
tering from long, thin dielectric cylinders. 
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Backscatter Coefficient of Smooth Surface 

The backscatter coefficient is determined by deriving 
an expression for the scattered electric field at the re- 
ceiver due to a plane wave incident on the surface. The 
scattered field is given in general by the Stratton-Chu 
integral as modified by Silver [21], 

The scattering geometry is illustrated in Figure IV-1. 
The unit vectors n^ and TI 2 define the incident and scat- 
tering-directions, respectively. Following Fung [22] and 
Leader [23] a local orthonormal coordinate system t, d 
is constructed as follows: 

d = 7i, 

The normal vector n is constant for the special case of a 
smooth surface. Furthermore, for the smooth surface 
t = X, where x is the unit vector along the x-axis. 

The general expression for the scattered field at a 
point P in the direction n^ is [23] 



//[(■>! xE) (^'‘H)]' 

ex p ( j r • ) c(S 
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where k * wave number 

R - distance of P from origin 

= intrinsic impedance of atmosphere 
r = position vector of surface element dS 

The cross products (ii x E) and (n x H) denote the local 
surface currents. Assuming that the incident field is 

= 5E„ exp(-jkv^, r) 

with time variation exp(jo)t) understood, the surface cur- 
rents may be expressed as 

(■n.-nKl -RJid'd) 

= +Rv)(a*a)(nxt) + 

(fi,-n){l t][-(-^)] 

where and R^ are the Fresnel reflection coefficients for 
horizontal and vertical polarizations, respectively. For 
horizontal polarization 

a = X 

and for vertical polarization 

{X - {j + Z 


t. E„ 


CIV-3) 



Therefore, the induced currents on a smooth surface are 



{««£)„ 

II 





= x{\ 


coxi 9 

for a 

horizontally polar iz< 

sd incident field 


{■n » E)v 

= -x(l 

-RvlE. 

coxiB 


(n * H)v 

= v(i 

*Rv)t 


for a 

vertically 

polarized 

incident 

field . 

these 

expressions 

into (IV 

-2) leads 

to 


Ej(p) = KEoi/'Ae-Kp[jkF-(^2-'nj]olS 

where 

^ _ -j kexp(-jkRj 

-^ttR 


CIV-4a) 


and 


CIV-4b) 


Substituting 


(IV-5) 


(IV-6) 


and A is a function of incidence and scattering angles and 
polarization. 
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A= T^j[-(n^*i})(l + Rj] 5(1 + Ry) 

CIV- 7a) 

- (fijX x)c^^(l -Ry) 

for vertical polarization, and 

A= X (l-Rjco<t^ - (n/ fi)(l + Rj] civ- 7 b) 

for horizontal polarization. 

For the special case of backscatter, 

yi^ ~ “H, = ’^IjAU/nO Z CotZ^ 

Thus, (IV- 5) becomes 

EjIp) = KEjA ff j kv- -(-2 ijMmO + 

2 7[ Cycni 0)]dx dij 

It is convenient to express the position vector as 

r = ~{r ^4A>n<k) X + ( r cotf <t>) y 

and to transform the variables of integration to polar 
coordinates. Thus, 
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Es(p) = KE,A/y%xp[jk(-2x»m9)' 

° (IV- 8) 

r 4>^r dr d4> 

Noting that Bessel functions of integer order are defined 
by [24] 

■ "It 

= 'hirijj. exp(jxcorf^)exp(j-n^)c(^ 

(IV- 8) may be written 

E^(p) = KEgA^ ' 2irJj-2krxt<^^) r dr 

(IV- 9) 

= KE/^'' 2 H( 2 k rMmO]rdr 


since J (x) is an even function, 
o 

Applying the theorem [24] 

fj dx = r'^J^ir] 

leads to 

F* / \ TrKEflA . . 


(IV- 103 


The radar cross section of the illuminated surface is 
determined by inserting (IV- 10) into (IV- 1). 
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a = 


E. 


I A I ^ 

L8k(W0)R 


CIV-11) 


2 

The area o£ illuminated surface is assumed to be . 

Therefore, the backscatter coefficient for the smooth sur- 


face is 




For horizontal polarization 


= 2R„ ood^ 


and for vertical polarization 


IaI = V(2FLc<«'S)^ (K^2e) 


Attenuation by Canopy 


(IV-12) 


That portion of the incident electromagnetic field 
which eventually reaches the soil experiences attenuation 
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as it propagates through the canopy. The attenuation con- 
stant is defined by (III-17). In terms of the field above 
the canopy, the electric field incident on the soil is 

r- r- -aH,<lec© 

E -, = Eo e 

where H is the height of the canopy and 0 is the angle of 
incidence . 

The effective backscatter coefficient of the soil is 
determined from the attenuation of power due to propaga- 
tion through the canopy. Since the incident electric field 
is reduced by exp(-aHsece) before reaching the soil, the 
incident power is reduced by exp (- 2aHsec0) . Furthermore, 
since the two-way path of the backscattered power is 
2Hsec0, the total attenuation of power due to the canopy is 
exp (- 4aHsec0) , Therefore, the effective backscatter coef- 
ficient of vegetated soil is 

/ o\ o riv- 1 ■^'1 

la = a e 

where o° is the backscatter coefficient of bare soil. 
Scattering by Canopy 

Two models for the scattering coefficient of a plant 
canopy have been considered. One is Peake's model for 
long, thin dielectric cylinders, discussed in Chapter II. 
The other is the Lommel-Seeliger model [25], which treats 
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the canopy as an ensemble of randomly oriented flat plates. 
Obviously, each model is intended to apply only to certain 
types of vegetation. 

The difficulty of applying either of these two models 
to the present development is that each one assumes the 
canopy has infinite depth. In other words, the scattering 
coefficient given by each model is the total average scat- 
tering coefficient of vegetated terrain. What is desired 
is a method of isolating the component of scattered radia- 
tion due to the canopy from the component due to the soil. 

A modification of Peake*s model seems to provide an 
appropriate solution to the problem. Peake's model is 
compatible with the model presently being developed because 
it accounts for exponential attenuation of the field within 
the canopy. It seems reasonable to assume that if the at- 
tenuation of the incident field is less than exp( -aHsec0) , 
the contribution of scattering from the soil is signifi- 
cant. As the canopy increases in height or density, scat- 
tering from the canopy becomes more significant and scat- 
tering from the soil becomes less significant. Such 
effects are accounted for by assuming that the component of 
scattering due to the canopy is 

o of 

q = cr^(/~e j 

where is the scattering coefficient computed from 
Peake's model. 



Theoretical Predictions 


For reference, the backscatter coefficient of a 
smooth soil surface is illustrated in Figure IV- 2. The 
backscatter coefficients for horizontal and vertical polar- 
izations are plotted as functions of incidence angle for 
two values of soil moisture content. An increase in mois- 
ture causes an increase in a®. For angles of incidence 
less than 25°, backscatter is about the same for both 
polarizations. For larger angles of incidence, backscatter 
is less for vertical polarization, but is more sensitive 
to variations of moisture. 

Assuming that the power scattered from the soil and 
from the canopy sum incoherently, the total backscatter 
coefficient of the vegetated surface is the sum of (IV- 13) 
and (IV-14). 


O' 


t 



(IV- 15) 


-4aHx«c0 -2aHxWC0 

= O' e + Opll - e 


The backscatter coefficient of (IV- 15) is plotted in 
Figure IV- 3 as a function of incidence angle for two values 
of soil moisture. The canopy is assumed to be 35 cm deep, 
and the incident radiation has a frequency of 1.5 GHz. At 
smaller angles of incidence an increase in moisture of 20% 
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Figure IV- 2. Backscatter coefficient of a smooth 
surface. 
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IV-3. Backscatter coefficient of smooth, 

uniformly vegetated surface as function 
of incidence angle. 
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causes an increase in the scattering coefficient of about 
5 dB. At incidence angles of about SO® for vertical polar- 
ization and 70° for horizontal, scattering from the canopy 
becomes more significant than that from the soil, causing 
a general increase in o*. Furthermore, as scattering from 
the canopy increases, the variation of o' due to moisture 
diminishes . 

The effect of changing the frequency of radiation is 
illustrated in Figure IV-4, in which the backscatter coef- 
ficient for frequencies of 3, 5, and 7 GHz is plotted. At 
higher frequencies (shorter wavelengths) scattering from 
the canopy becomes predominant at smaller incidence angles. 
For a single frequency the same effect results from in- 
creasing the height of the canopy. 

Figure IV- 5 illustrates the dependence of the back- 
scatter coefficient on soil moisture content for three 
frequencies. The angle of incidence is 30°. At 1,5 GHz 
an increase in moisture of 20% causes an increase in the 
scattering coefficient of about 5 dB. At 3.0 GHz the 
change in scattering coefficient is 1.6 dB , and at 5.0 GHz 
the change is negligible. At 5.0 GHz scattering from the 
canopy is great enough to cause a general increase in a ' . 

Row Crops 

Since the contributions to scattering due to the 
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Figure IV-5. Dependence of o' on soil moisture 
content for various frequencies. 
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canopy and. to the soil have been determined in general, the 
determination of the average backscattering coefficient of 
row crops is primarily a geometrical problem. The results 
are similar to those obtained for the average apparent tem- 
perature of row crops in Chapter III. 

Backs catter Perpendicular to Rows 

For the special case of incidence and scattering per- 
pendicular to the rows, three scattering regions are de- 
fined as in Figure IV-6. In regions I and II both the soil 
and the canopy contribute to the scattered field, whereas 
in region III only the soil causes scattering. The angle 
of incidence in region I is (tt/ 2 - 6) rather than 6 as in 
regions II and III. 

The component of the backscatter coefficient due to 
the canopy is a weighted average of the scattering coeffi- 
cients in regions I and II, The area upon which the field 
is incident in region I is projected onto a plane parallel 
to the surface. Thus, the effective backscatter coeffi- 
cient of region I is 

cr;(«) = 

The total effective component due to tie canopy is 

-e)H * m’M i, 

^ B + W ' r (IV- 16 ) 
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Figure IV-6. Dimensions o£ rows and basic regions 
of backscatter. 
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where d is the average length of propagation through the 
canopy, as defined by (111-23). 

The expression of (IV- 16) is not applicable in gen- 
eral, since it was derived for the case in which B < A. 
However, it can be shown that generality is obtained by 
substituting the function L(A,B) , defined in Chapter III 
(p. 65), for B and L(A,B)cot6 for H. Thus, 


UA.BI w ' ' 


7) 


The total backs cattering coefficient is the average 
of all scattering components in regions I, II, and III. 


a = 


_i 

A+W 


[UA.Bl-^wKor'.aV'”") 


(IV-18) 


-h 


[A-L|A.B)]a' 


where a® is given by (IV-17) and o° is given by (IV-12). 
Backs cattering at Any Azimuth Angle 

Since the expression for scattering perpendicular to 
the rows is identical in form to the expression for appar- 
ent temperature, stated as (III- 30), it is reasonable to 
assume that the generalization for arbitrary azimuth an- 
gles, given in (III-28), is also applicable in the present 
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case. Therefore, it is stated without explanation that the 
only modification required is to define the distance B as 

B = H ’tam.O 


Theoretical Predictions 

The relation of the backscatter coefficient to soil 
moisture content depends on several factors. In general 
it may be said that as the component of scattering due to 
vegetation increases, the dependence on soil moisture de- 
creases, However, the effect of vegetation can be in- 
creased by changing frequency, canopy dimensions, or angle 
of incidence. 

The effect of frequency on the backscatter coefficient 
is illustrated in Figure IV- 7, in which three curves are 
plotted as functions of soil moisture. The plant row is 
assumed to be 50 cm high and SO cm wide. The angle of in- 
cidence is 30°, and the azimuth angle is 45°. The data 
indicate that at 1.0 GHz an increase in moisture of 20% 
causes the scattering coefficient to rise 4.0 dB, while at 
5,0 GHz the change is negligible. 

Figure IV-8 illustrates that the same effect is 
caused by changing the dimensions of the canopy. The data 
have been computed for a frequency of 3.0 GHz. For a 



Backscatter Coefficient (dB) 



Backscatter Coefficient (dB) 


104 



Figure IV-8. Dependence of o' on soil moisture content 
for various amounts of plant cover. 
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canopy height of 25 cm an increase in moisture of 201 
causes an increase in the backscattering coefficient of 
3.5 dB. For a canopy height of 75 cm the same increase in 
moisture causes a change in the scattering coefficient of 
on ly 1.5 dB . 

The angle of incidence partly determines the effect of 
vegetation on scattering, as illustrated in Figure IV-9. 
Each curve in this figure represents the backs catter coef- 
ficient, as a function of soil moisture, for incidence per- 
pendicular to the rows. At smaller angles of incidence, 
the dependence of the scattering coefficient on soil mois- 
ture is appreciable. However, at larger angles of inci- 
dence, in which case more of the soil is shadowed, the 
scattering coefficient becomes practically independent of 
soil moisture. 

Scattering from a Rough Surface 
Surface Model 

Johnson^s model of a rough surface, defined by a 
surface roughness factor, 6, and a coefficient of effective 
area, 3, is applicable to the problem of scattering from a 
vegetated rough surface. The derivation of the model, with 
application to thermal emission, is discussed in Chapter 
III. In the present section only those points which per- 
tain especially to the scattering problem are discussed. 
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Figure IV- 9. Dependence of a' on soil moisture 

content for various angles of incidence. 


According to Johnson’s derivation [7], the incident 


and scattered radiation are related as follows: 




00<l9i 

| 3 hc^ 9 ^^ 

C0<i9i 




(IV- 19) 


where and are the scattered power densities for ver- 

tical and horizontal polarizations, ly and are the inci 
dent power densities, 6^ is the angle of incidence, and 
is the angle between the incidence direction and the 
effective normal of the scattering surface element* The 
scattering coefficients are defined as 

'^vv ~ Pv Ph 

7„h “ Ph ^ Pv 

(IV-20) 

Co^CLf - cxTtf a- ✓(tvno.j. 

T'vh ~ P^ o-j - p^ycu^ci. 

where and are the Fresnel reflection coefficients 

and and are the angles of rotation of incident and 
reflected polarization fields into the local plane of 
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incidence. 

Since only backscattering is being considered, the 
local angle of incidence must be zero. In other words, the 
effective normal of the scattering surface element must be 
parallel to the direction of propagation of the incident 
radiation. Consequently, 0^^ = 0 for all incidence angles. 
Furthermore, since backscatter is assumed to be due only to 
specular reflection of normally incident fields, there is 
no depolarization of the incident fields. Therefore, the 
scattering coefficients given by (IV-20) reduce to 



(IV-21) 




From (IV- 19) the backscatter coefficients of the rough soil 
are 



Ph 2 
oo<i 9 


(IV-22a) 


for vertical polarization and 


0 


9 


(IV-22b) 


for horizontal polarization. 
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Theoretical Predictions 


The backscatter coefficient of row crops on a rough 
surface is given by [1^-18) with defined by (IV-22]. 

To evaluate the surface roughness factor, 6, and the 
coefficient of effective area, 3, the surface must be de- 
scribed by the distribution of the surface normal. The 
derivation of Johnson's model assumes that the normal is 
distributed uniformly with respect to the azimuth angle. 
Therefore, the specified probability density function is a 
function only of the angle of incidence. A suitable func- 
tion for this application is 

f {^) = 

which leads to 

6 - 1.178 


and 


2 ootf 


The backscatter coefficient is plotted as a function 
of incidence angle in Figure IV-10 for two values of soil 
moisture. A 5 GHz signal is assumed to be incident per- 
pendicular to plant rows 60 cm high and 40 cm wide. It is 
clear that a 20% increase in moisture causes an increase in 
o' of about 4 dB for both horizontal and vertical 



Backscatter Coefficient (dB) 


no 



Figure IV-10, Backscatter coefficient of a row crop 
on a rough surface. 
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polarizations at small angles o£ incidence. However, at 
larger angles of incidence scattering is less sensitive to 
variations in moisture. 

The effect of varying the azimuth angle is illustrated 
in Figure IV-11. There is no noticeable change in the 
moisture sensitivity of the scattering coefficient for 
horizontal polarization. However, for vertical polariza- 
tion, the sensitivity to moisture is greater at large 
azimuth angles (approximately parallel to rows). 

The sensitivity of the backscatter coefficient to 
variations of moisture depends upon the operating frequen- 
cy, as illustrated in Figure IV-12. At 1.0 GHz an increase 
in soil moisture of 20^ causes an increase in the scat- 
tering coefficient of 5 dB. At 5.0 GHz the change in the 
scattering coefficient is 3.5 dB, and at 9.0 GHz it is only 
1.4 dB. 

The same effect results from changing the dimensions 
of the plant rows, as illustrated in Figure IV- 13. A 5 GHz 
signal is assumed to be incidence on plant rows of three 
different sizes. The backscatter coefficient of 38 x 33 
cm rows increases 3.83 dB with an increase in soil moisture 
of 20%. For 100 x 70 cm rows the change in scattering 
coefficient is 1.12 dB. 

As expected, for any particular frequency or vegeta- 
tion condition, the backscatter coefficient of a rough 



Backscatter Coefficient (dB) 
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surface is more sensitive to soil moisture than that of a 
smooth surface, since the contribution of backs cattering 
from the soil is greater. 
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CHAPTER V 

COMPARISON OF THEORETICAL PREDICTIONS 
WITH EXPERIMENTAL DATA 

Measured apparent temperatures and backscatter 
coefficients of vegetated terrain are presented in this 
chapter to illustrate some of the effects predicted by the 
models developed in Chapters III and IV. Unfortunately, 
none of the data, much of which is taken from the litera- 
ture, -is accompanied by documentation of both soil moisture 
and vegetation. Therefore, apparent temperatures and back- 
scatter coefficients cannot be illustrated as functions of 
soil moisture. However, the effects of other parameters 
on scattering and emission can be illustrated. 

A program is presently being conducted at Texas A^M 
University to measure the apparent temperatures of test 
sites for which soil moisture, surface roughness, and vege- 
tation are well documented. However, no data are yet 
available . 

The increase of apparent temperature due to vegetation 
is clearly illustrated in Figure V-1. The fields pictured 
lie along the flight line near Weslaco, Texas. Data from 
the March 14, 1972, flight over these fields were provided 
by NASA/Goddard Space Flight Center. The data plotted are 
from a 1.42 GHz dual-polarization radiometer and a 19.4 



APPARENT TEMPERATURE 



WESLACO, TEXAS march 1972 



Figure V-1. Apparent temperatures of bare and vegetated fields at 
Weslaco, Texas . 
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GHz scanning radiometer. The center line o£ the digitized 
scanner image is plotted. The apparent temperature is as 
much as 65 degrees higher in vegetated fields than in bare 
fields. This increase is apparently due only to the vege- 
tation, since all fields are believed to have had approxi- 
mately equal moisture levels. 

A change in dependence on polarization, as predicted 
by Figure III-12 (p. 67), is also apparent in Figure V-1 
Cp* 117). Although the vertical channel of the 1.42 GHz 
data is consistently higher than the horizontal channel, 
the difference between the two channels is considerably 
less for vegetated fields than for bare fields. 

The lack of dependence on polarization is better 
illustrated by Figure V-2, in which apparent temperatures 
of alfalfa and oats are plotted as functions of refraction 
angle. These data and other data presented below were 
obtained by Peake and Oliver [12]. The apparent tempera- 
tures were measured at a frequency of 10 GHz. It is clear 
that for any angle, a change of polarization causes a 
change in apparent temperature of no more than about 5°C. 
The peak at the Brewster angle, characteristic of vertical 
polarization, is suppressed. 

The difference between apparent temperatures of wet 
and dry fields is illustrated in Figure V-3. The apparent 
temperatures of two soybean fields, one of which was irri- 
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Figure V-3. Apparent temperature of green soybeans at X-band [12] 
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gated, are plotted as functions of refraction angle. It 
seems that the apparent temperature of the wet field is 
significantly lower than that of the dry field for angles 
greater than 10°, However, it must be noted that the 
actual surface temperatures differ by 7°C. Raising the 
solid curves to compensate for this difference, it is seen 
that the difference in apparent temperatures is actually 
small at larger refraction angles. There seems to be 
greater sensitivity to soil moisture at the lower angles , 
as predicted in Chapter III. 

The backscatter coefficient of vegetated terrain 
is almost independent of polarization (assuming trans- 
mitting and receiving polarizations are the same), as 
illustrated in Figure V-4. In this figure the backscatter 
coefficient of green oats, measured at X-band is plotted 
as a function of incidence angle. A change of polariza- 
tion causes a change in the scattering coefficient of no 
more than 2 . 5 dB . 

The frequency dependence of the backscatter coeffi- 
cient is illustrated in Figure V-5. The data plotted are 
backscatter coefficients of green oats, 25 cm in height, 
measured at frequencies of 1.8, 10, and 35 GHz, The 1.8 
GHz data are from 3 to 6 dB lower than the higher fre- 
quency data, whereas the data for 10 and 35 GHz differ by 
no more than 1 dB . These data imply that at 1.8 GHz a 
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Figure V- 5. Backscatter coefficient of green oats 
25 cm high at various frequencies. 


significant part of the backscatter is due to the soil, 
whereas at frequencies higher than 10 GHz backscatter is 
due almost entirely to the vegetation. 

The fact that increasing canopy height has the same 
effect as increasing frequency is illustrated in Figure 
V-6. Backscattering coefficients of three fields of wheat 
9 cm, 36 cm, and 74 cm in height, are plotted as functions 
of incidence angle for frequencies of 5.87 GHz and 9.375 
GHz. At the lower frequency, data for the 74 cm canopy 
show little variation with incidence angle. In contrast, 
the data for the two shorter canopies have steep slopes 
between 0® and 10®, characteristic of a smooth surface. 

At the higher frequency, the data for the 9 cm canopy stil 
exhibit a sharp decrease between 0® and 10® incidence. 
However, the data for the 36 cm canopy tend to look more 
like the data from the tallest canopy. In effect, the 
higher frequency makes the canopy appear more heavily 
vegetated , 

From the data presented in this chapter it is not 
possible to assess the quantitative accuracy of the models 
developed in Chapters III and IV. However, the results of 
these models have been shown to be at least qualitatively 
realis tic . 
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CHAPTER VI 
CONCLUSION 

In Chapters III and IV, models are developed for the 
apparent temperature and backscatter coefficient of vege- 
tated terrain. Data computed from these models indicate 
that the sensitivity of the apparent temperature and back- 
scatter coefficient to variations of soil moisture depends 
not only on the amount of vegetation, but also on certain 
system parameters. The various factors which help deter- 
mine the sensitivity of data to soil moisture are 


(1) 

canopy height 


(2) 

canopy density 


(3) 

row width (for row 

crops) 

(4) 

row spacing (for row crops) 

(5) 

signal frequency 


(6) 

angle of incidence 


(7) 

azimuth angle (for 

row crops) 


Since several parameters have similar effects on the 
data, it should be expected that numerous equivalent states 
exist. Equivalent states are defined as different combina- 
tions of the above parameters which cause equivalent appa- 
rent temperatures or backscatter coefficients. By holding 
the system parameters constant, the number of variables is 
reduced to two (for 100% coverage) or four (for rows) . It 
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is desirable to find a single parameter by which equivalent 
vegetation states may be identified. 

Several equivalent states, determined by . inspection 
of calculated apparent temperatures, are listed in Tables 
VI - 1 through VI-3, The data represented by the tables are 
for frequencies of 1, 3, and 5 GHz. In each table the 
vegetation states are defined by canopy height and vege- 
tation density, A uniform canopy (1001 coverage) is 
assumed. Vegetation density is defined as the volume of 
plant matter per unit volume of space containing the can- 
opy. Comparison of the tables indicates that equivalence 
of states is dependent on frequency. However, for a uni- 
form canopy, equivalence is independent of the angle from 
which the terrain is observed. 

Also listed in Tables VI-1 through VI- 3 are the pro- 
ducts of height and density of each state. In Table VI - 1 
it is seen that the equivalent states having densities of 
0.04 and 0.05 have approximately equal height-density pro- 
ducts. It is seen in Tables VI-2 and VI - 3 that all equiva- 
lent states have approximately equal height-density pro- 
ducts., Therefore, the height- density product may be used 
to identify any of the equivalent states. 

The model developed in Chapter III predicts that the 
presence of vegetation decreases the sensitivity of the 
apparent temperature to variations of soil moisture by 
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TABLE VI -1 

EQUIVALENT VEGETATION STATES AT 1 GHZ 


Vegetation Density 


0.03 

0 

.04 

0 

1.05 

Canopy 

Height 

(cm) 

Height- 

Density 

Product 

Canopy 

Height 

(cm) 

Height- 
Density 1 
Product 

Can opy 
Height 
(cm) 

Height- 

Density 

Product 

50 

1.50 

29 

1.16 

23 

1.15 

60 

1.80 

35 

1. 40 

28 

1.40 

70 

2. 10 

41 

1.64 

33 

1.65 

80 

2. 40 

46 

1. 84 

38 

1,90 

90 

2. 70 

51 

2.04 

42 

2.10 

100 

3.00 

56 

2.24 

47 

2. 35 



70 

2. 80 

60 

3.00 



80 

3.20 

68 

3.40 



90 

3.60 

77 

3. 85 
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TABLE VI -2 

EQUIVALENT VEGETATION STATES AT 3 GHZ 


Ve ge t at i on Dens i ty 


0.03 

0.04 

0.05 

If 

Height- 

Density 

Product 

Canopy 

Height 

(cm) 

Height- 

Density 

Product 

Canopy 

Height 

(cm) 

Height- 
Dens ity 
Product 

40 

1.20 

30 

1.20 

24 

1.20 

50 

1.50 

38 

1.52 

31 

1.55 

60 

1.80 

45 

1.80 

37 

1. 85 

70 

2.10 

53 

2.12 

43 

2.15 

80 

■ 2. 40 

61 

2.44 

49 

2.45 

90 

2. 70 

68 

2. 72 

55 

2. 75 

100 

3.00 

76 

3.04 

62 

3.10 



90 

3.60 

73 

3.65 



100 

4,00 

82 

4. 10 
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TABLE VI -3 

EQUIVALENT VEGETATION STATES AT 5 GHZ 


Ve ge t ati on De ns i ty 


0.03 

0.04 

0.05 


Height- 

Density 

Product 

Canopy 

Height 

(cm) 

Height- 

Density 

Product 

Canopy 

Height 

(cm) 

Height- 

Density 

Product 

50 

1.50 

40 

1.60 

32 

1.60 

60 

1.80 

46 

1.84 

38 

1.90 

70 

2.10 

54 

2.16 

44 

2.20 

80 

2. 40 

62 

2.48 

49 

2.45 

90 

2. 70 

69 

2.76 

56 

2. 80 

100 

3,00 

76 

3.04 

63 

3.15 



90 

3.60 

73 

3.65 



100 

4.00 

83 

4.15 
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increasing the apparent temperature with respect to that of 
bare soil. Figure VI - 1 illustrates the difference between 
the apparent temperatures of bare and vegetated soil for 
various frequencies and vegetation states. All data are 
computed for 0=0. A soil moisture content of 20% is as- 
sumed. As expected, the difference in apparent tempera- 
tures is relatively small at 1 GHz, whereas it is quite 
significant at higher frequencies. For a height-density 
product of 4, the temperature difference is 5°K at 1 GHz, 
39°K at 3 GHz and 70°K at 5 GHz. 

The effect of differences in apparent temperatures 
of bare and vegetated soil on the ability to determine 
moisture content is illustrated in Figure VI-2. The mois- 
ture content plotted is that which would cause bare soil to 
have the same apparent temperature as vegetated soil with a 
moisture content of 20%. It is clear from this figure that 
vegetation causes the soil to appear to have a lower mois- 
ture content. For a vegetation height- dens ity product of 
2, the soil appears to have a moisture content of 19% at 
1 GHz, 14% at 3 GHz, and 4% at 5 GHz. 

Inspection of backscatter coefficients computed from 
the model developed in Chapter IV reveals that the height- 
density product of vegetation is a good indicator of equiv- 
alent states for backscatter. However, for a smooth sur- 
face, equivalence depends on angle of incidence. It is 




Measured Moisture Content (Percentj 



Figure VI-2. Moisture content determined from 

apparent temperature of vegetated soil 
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found that vegetation states having the same height-density 
product have backscatter coefficients differing by no more 
than 10% for incidence angles less than 30®. For such 
angles, equivalence may be assumed when height-density pro- 
ducts are equal. However, for angles of incidence greater 
than 30°, equal height- density products do not necessarily 
imply equivalence of states. 

Plant rows are more difficult to describe than a uni- 
form canopy. In addition to the parameters already con- 
sidered, row width, row spacing, and azimuth angle must be 
considered. The height-density product is not adequate to 
identify a vegetation state, even if azimuth angle and row 
spacing are held constant. Nevertheless it is expected 
that a suitable factor can be found to identify equivalent 
states of row vegetation. Having determined such an indi- 
cator, curves such as those plotted in Figures VI-1 (p. 132 ) 
and VI-2 (p . 133 ) can be plotted for a number of other 
vegetation states. 

Plots of measured moisture content vs. vegetation 
state make it possible to improve the accuracy of remote 
soil moisture measurements of vegetated terrain. However, 
such curves are difficult to produce in practice, since 
many parameters must be known. Obviously the vegetation 
state must be identified by a suitable indicator (such as 
the height-density product). Furthermore, the actual 
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surface temperature must be known. For all calculations 
in this report the surface temperature is assumed to be 
298®K. In addition, surface roughness and atmospheric 
effects must be accounted for. 

There is clearly a need for more experimental data 
from terrain for which soil moisture and vegetation are 
known, little of which is available at present. The re- 
sults of this investigation indicate that the effect of 
vegetation on microwave data from terrain is greater than 
has been anticipated in some previous work. From the theo- 
retical and experimental data presented here, it is clear 
that the ability to measure soil moisture by remote sensing 
techniques is highly dependent on the amount of vegetation 
over the soil. However, the sensitivity of data to varia- 
tions of soil moisture can be improved by adjusting certain 
system parameters. For example, airborne radiometer data 
have already verified that low frequency sensors are more 
sensitive to variations of soil moisture [11]. Ground- 
based experiments are needed to determine critical amounts 
of vegetation, above which soil moisture cannot be deter- 
mined, for various frequencies. If data from future air- 
borne experiments are to be of any value for soil moisture 
measurement, efforts must be made to compensate for vege- 
tation. It is hoped that the result of such efforts will 
be to improve the techniques by which soil moisture content 
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is measured by microwave remote sensors. 
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APPENDIX A 

PENETRATION OF MICROWAVES INTO SOIL 

The depth of penetration, or skin depth, of a medium 
is defined as the distance over which an electric field is 
attenuated to 1/e (37%) of its initial strength as it 
propagates through the medium [17]. Since the intensity of 
a field propagating in the x-direction is 

E(x)=E,e'“’^ 

it is easy to see that the skin depth is equal to the reci- 
procal of the attenuation constant. In terms of the elec- 
trical properties of the medium, the skin depth is 

S = 

where to = radian frequency of signal 

p = magnetic permeability of medium 
€ = dielectric constant of medium 
a = conductivity of medium 

The skin depth of soil is a function of the soil mois- 
ture content, as illustrated in Figure A-1. In this figure 
the skin depth of three types of soil, the permittivity of 
which was measured by Wiebe [7], is plotted as a function 
of soil moisture for a frequency of 9 GHz. For these three 




141 


Gila Sandy Loam 

Tarrant Stony Clay 



Figure A-1. Skin depth of various soils as a 
function of moisture content. 


Skin Depth (Wavelengths) 
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soils the skin depth is greater than half of a wavelength 
at moisture contents less than 10^, and less than half of 
a wavelength at moisture contents greater than 151. 

From (A-1) it is clear that the skin depth is also a 
function of frequency. The skin depth of Gila Sandy Loam 
is plotted in Figure A- 2 for frequencies of 1,5 GHz and 
9.0 GHz. The permittivity of the soil is assumed to be 
the same at both frequencies. At the lower frequency the 
skin depth is greater than 4 cm for all values of moisture 
content up to 201, At the higher frequency the skin depth 
is less than 2 cm for any moisture content greater than 
10 %. 

A knowledge of the skin depth of soil at various fre- 
quencies is helpful in assessing the value of passive 
microwave sensors for remote measurement of soil moisture. 
The skin depth is an indication of the depth from which 
radiation contributes significantly to the apparent tem- 
perature. Since the surface of the ground is usually dry, 
regardless of subsurface moisture content, it is desirable 
that the apparent temperature include a significant con- 
tribution from several centimeters beneath the surface. 
Therefore, longer wavelength sensors are advisable for 
this application. 



Skin Depth (cm) 
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APPENDIX B 

ATMOSPHERIC EFFECTS ON APPARENT 
SURFACE TEMPERATURE 

The atmosphere acts in several ways to alter the 
electromagnetic signal which radiates from the ground. 

Some of the diffuse radiation from the atmosphere is scat- 
tered by the target into the antenna beam, augmenting the 
signal due to thermal radiation from soil and vegetation. 
As the combined radiated and scattered energy propagates 
toward the antenna, it is attenuated by the intervening 
atmosphere. At the same time, however, the atmosphere 
between the target and the sensor is radiating energy 
directly into the antenna, adding another component to 
the total apparent temperature. 

At microwave frequencies the primary causes of atmos- 
pheric attenuation are water vapor and oxygen. The atten- 
uation due to each of these gases is illustrated in 
Figure B-1. These curves are computed using a modified 
form of Van Vleck^s model [26]. The attenuation curve 
for water vapor has maxima at 22.3 GHz and about 180 GHz. 
(The higher frequency is not included in the figure.) 

The curve of oxygen absorption has maxima at 60 GHz and 
120 GHz. 

Summing the attenuation due to water vapor and oxygen 



Attenuation (dB/kiri) 
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along with a third term to account for absorption by other 
gases, gives the total atmospheric attenuation, as illus- 
trated in Figure B-2. The peaks at 20 and 60 GHz are still 
apparent. It is clear that atmospheric attenuation de- 
creases constantly below 20 GHz, becoming negligible below 
about 1 GHz, 

The actual effect of atmospheric attenuation and 
radiation is better understood by evaluating the apparent 
temperature of a theoretical rough surface at different 
frequencies. Figure B-3 illustrates apparent temperatures 
computed from Peake's model as modified by Richerson [27]. 
The sensor is assumed to be at an altitude of 609.6 m 
(2000 ft) above the surface. Data were computed for 
various values of humidity at frequencies of 1.42, 10.69, 
19,35, and 31,4 GHz. It is clear that the effect of 
changing humidity is negligible at 1.42 GHz, but becomes 
more pronounced at higher frequencies. The number at the 
right-hand end of each curve is the apparent temperature 
which would be measured if the atmosphere only attenuated 
the signal, radiating no energy itself. At the lower 
frequencies the atmosphere adds about 2 degrees to the 
measured apparent temperature, whereas at higher fre- 
quencies the contribution of the atmosphere may be 9 or 
10 degrees. 




Apparent Temperature (°K) 
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Figure B-3. Atmospheric effects on measured apparent 
temperature at various frequencies. 



